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I.  INTRODUCTION 


This  paper  presents  the  results  of  Monte  Carlo  studies  that 
have  been  performed  to  determine  if  the  presence  of  dust  within  a cloud 
at  altitudes  between  14  and  15  km  could  be  detected  at  satellite  alti- 
tudes. Although  the  purpose  of  the  present  study  is  to  determine  what 
wavelength  intervals  would  be  best  for  detection  purposes,  all  of  the 
Monte  Carlo  studies  that  have  been  performed  to  date  have  been  in  the 
2.7  y®  wavelength  band.  In  addition,  some  calculations  have  been  run 
with  the  MIE-2  procedure  to  obtain  scattering,  absorption,  and  extinction 
cross  sections  for  both  the  cloud  and  the  dust  for  wavelengths  between 
2 and  15  pm  to  assist  in  determining  what  wavelength  Intervals  would  be 
most  promising  for  use  in  detection  of  dust-loaded  clouds. 

This  work  was  perfomed  under  Contract  Numbers  06-0481  and 
06-6328  with  the  Sandia  Laboratories.  Dr.  R.  D.  Rawcliffe  of  the  Aero- 
space Corporation  has  contributed  to  the  study  by  both  suggesting  the 
study  and  by  defining  the  parameters  to  be  used  to  describe  both  the 
cloud  and  the  dust  to  be  loaded  to  the  cloud. 
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II.  CALCULAT IONAL  METHOD 


The  basic  computational  tools  being  used  in  the  study  are  the 
MIE-2  procedure  (Ref.  1)  and  the  FLASH  Monte  Carlo  procedure  (Ref.  2). 

The  MIE-2  procedure  uses  MIE  theory  for  spherical  particles  with  a 
complex  index  of  refraction  and  a defined  particle  size  distribution 
to  compute  phase  function  data  and  absorption,  scattering,  and  extinction 
cross  section  data. 

FLASH  is  a Monte  Carlo  program  which  treats  the  scattering  and 
absorption  of  both  sunlight  and  air  and  ground  emission  in  a spherical- 
shell  atmosphere.  FLASH  is  capable  of  tieating  temperature  and  pressure 
dependent  absorption  in  the  infrared  by  C02,  H20,  C>3,  N20,  CO,  CH^,  and 
02  using  the  gaseous  absorption  cross  section  data  generated  by  McClutchey 
and  Shelby  at  the  Air  Force  Geophysical  Laboratory  (Ref.  3).  The  restric- 
tion to  these  seven  gases  in  FLASH  is  due  only  to  the  unavailability  of 
data  for  other  gases  in  the  AFGL  data  base.  The  AFGL  data  are  averaged 
over  a wavenumber  interval  of  20  cm  * and  are  available  at  5 cm  * spacings 
for  wavenumbers  between  10  cm  and  13,155  cm-1  (0.76  to  1000  pm) . The 
AFGL  data  are  input  to  FLASH  for  seven  temperatures  between  175°  and 
325°K.  Data  for  temperatures  between  those  used  in  the  data  base  are 
obtained  in  FLASH  through  Interpolation. 

When  using  the  AFGL  data  in  FLASH,  it  is  required  that  the  wave- 
number  range  of  interest  be  divided  into  a number  of  intervals  that  are 
200  cm  ^ or  less  in  width  over  which  the  scattering  and  continuum 
absorption  cross  sections  can  be  assumed  to  not  vary  with  wavenumber  in 
each  interval.  Within  each  of  these  200  cm  ^ intervals,  the  effects 
of  gaseous  absorption  and  scattering  can  be  calculated  for  up  to  10 
sub intervals  of  up  to  20  cm  width. 

The  Initial  effort  during  the  first  three  months  of  the  study 
was  Involved  in  making  the  FLASH  procedure  operational  on  the  Sandia 
CDC-6600  computer  system.  FLASH  problems  were  then  run  for  a model  mid- 
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latitude  summer  atmosphere  containing  a 15  mean-free-path  thick  cumulus 
cloud  located  between  14  and  15  km  altitude  and  loaded  with  different 
amounts  of  dust.  The  wavelength  range  considered  in  the  calculations 
uses  2.653  pm  X 2.760  pm. 

The  atmospheric  model  used  for  these  calculations  was  for  a 
25-km  meteorological  range  in  a mid-latitude  summer  atmosphere.  The 
variation  of  the  water  vapor  content  and  the  temperature  and  pressure 
with  altitude  was  obtained  from  measurements  taken  on  27  July  1973 
that  were  previously  supplied  to  RRA  by  the  Air  Force  Technical  Appli- 
cations Center  under  another  project.  A 15  mean-free-path  thick  cumulus 
cloud  for  2.7  pm  wavelength  light  was  positioned  at  altitudes  between 

14  and  15  km. 

The  altitude  dependence  of  the  aumospheric  temperature  as  used 
in  the  FLASH  calculations  is  shown  in  Fig.  1.  It  is  seen  in  Fig.  1 
that  the  atmospheric  temperature  has  a minimum  (~205°K)  between  14  to 

15  km  altitude.  It  should  also  be  noted  that  above  the  cloud  the 
temperature  reaches  a maximum  at  an  altitude  of  50  km. 

The  particle  size  distribution  used  for  the  dust  in  the  cloud 
was  that  described  by  the  equation 

^ - 1.58468r'4,5  (p"1  cm-3) 
dr 

for  Sahara  dust  (Ref.  4).  The  value  of  the  constant  in  the  size  distri- 

3 

butlon  was  determined  by  assuming  that  the  volume  of  dust  in  a one-m 
-5  3 3 

volume  was  5X10  cm  /m  . This  resulted  in  a value  of  N ■ 8569.15 
3 

particles/cm  for  dust  particles  with  radii  less  than  10  p.  The  index 
of  refraction  at  a wavelength  of  2.7  pm  was  taken  to  be  1.46  - 0.0251 
(Ref.  4). 

The  cumulus  cloud  was  assumed  to  have  a particle  size  distribu- 
tion of  the  form  (Ref.  5) 
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The  cloud  droplets  were  assumed  to  be  water  and  to  have  an  index  of 

refraction  of  1.134  - 0.005049i  (Ref.  6)  for  a wavelength  of  2.7  um. 

3 

N was  taken  to  be  100  particles/cm  . The  parameter  a was  determined 
by  integrating  the  size  distribution  over  the  particle  radius.  The 
optical  thickness  at  2.7  um  of  a one-km  thick  cumulus  cloud  was  then 
determined  from  MIE  calculations  to  be  22.166.  This  value  was  multi- 
plied by  0.676712  to  give  an  optical  thickness  of  15  for  the  cloud. 

The  extinction,  scattering,  and  absorption  optical  thickness 
for  the  dust-loaded  cloud  are  shown  in  Table  I.  These  optical  thick- 
nesses do  not  include  the  absorption  optical  thickness  resulting  from 
gaseous  absorption  in  the  one-km  thickness  of  the  cloud.  The  optical 
thickness  listed  in  Table  I are  for  1)  cloud  with  no  dust,  2)  cloud 
with  standard  dust,  3)  cloud  with  100  times  standard  dust,  and  4) 
cloud  with  1000  times  standard  dust,  where  standard  dust  is  defined  by 
the  Sahara  dust  distribution  described  above.  Adding  the  standard  dust 
content  to  the  cloud  resulted  in  a negligible  change  in  the  extinction, 
scattering,  and  absorption  optical  thicknesses  of  the  cloud  region. 

Adding  100  times  and  1000  times  standard  dust  caused  noticeable  increases 
in  the  optical  thicknesses  of  the  cloud  region. 

The  geometry  model  used  in  the  FLASH  calculations  is  shown  in 
Fig.  2.  The  receiver  is  positioned  at  satellite  altitude  and  the  posi- 
tion viewed  on  the  top  of  the  cloud  is  defined  in  terms  of  the  zenith 
angle-of-view,  0,  and  the  azimuthal  angle  <t>  between  the  plane  contain- 
ing the  line-of-sight  and  a radial  to  the  satellite  and  a plane  con- 
taining the  radial  to  the  satellite  and  the  radial  to  the  sun.  The 
angle  0^  is  the  zenith  direction  to  the  sun  at  the  sub-satellite  point 
and  the  angle  CXq  is  the  solar  zenith  angle  at  the  point  of  observation 
at  the  top  of  the  cloud. 

The  FLASH  calculations  for  incident  solar  radiation  were  run 
for  solar  zenith  angles  of  0^  “0°,  10°,  30°,  50°,  and  80°.  The  satellite 
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Fig.  2.  Geometry  Used  in  Monte  Carlo  Calculations 


altitude  was  taken  to  be  35,800  km.  Receiver  zenith  view  angles  of 
0 ■ 171.99,  173.77,  175.55,  177.33,  and  179.11°  and  azimuthal  angles 
of  $ ■ 0°,  90°,  and  180°  were  used  to  define  the  positions  viewed  by  the 
satellite. 

The  FLASH  procedure  is  a backward  Monte  Carlo  program  where 

photon  histories  are  started  at  the  receiver  position  in  each  of  the 

defined  directions  being  viewed.  The  results  are  obtained  in  units  of 
-2  -1  -2 

photons  cm  sr  per  photon  cm  in  a plane  normal  to  the  direction 

to  the  sun.  Although  the  results  in  this  paper  for  incident  sunlight 

will  be  shown  in  those  units,  the  results  can  be  converted  to  absolute 

power  units  for  a wavelength  range  interval  of  2.653  ym  £ X £ 2.76  urn 

by  multiplying  the  FLASH-generated  data  by  the  solar  constant  for  that 
-3  2 

interval  of  4.3X10  watts/cm  ym.  In  the  FLASH  calculations,  this 
wavelength  interval  was  represented  by  seven  wavenumber  intervals  of 
approximately  20  cm  ^ in  width.  For  the  sunlight  scattering  calcula- 
tions, the  incident  sunlight  was  assumed  to  be  uniformly  distributed 
in  the  total  wavelength  interval.  For  the  air  emission  calculations, 
the  source  strength  was  dependent  on  the  atmospheric  temperature  and  the 
absorption  cross  section  variation  with  altitude  for  each  of  seven  wave- 
number  intervals  of  approximately  20  cm  ^ in  width  in  the  wavelength 
range  between  2.653  ym  and  2.76  ym. 
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III.  CALCULATIONAL  RESULTS 


The  dependence  of  the  magnitude  of  the  scattered  solar  radiation 
at  the  satellite  on  the  number  of  collisions  followed  in  the  Monte  Carlo 
calculations  was  examined  to  insure  that  the  calculations  were  carried 
out  to  a large  enough  number  of  collisions  so  that  the  results  would 
not  be  underestimates  of  the  scattered  radiation  intensity.  Figure  3 
shows  results  typical  of  that  obtained  when  the  cloud  did  not  contain 
any  dust.  The  results  shown  are  for  a receiver  zenith  angle  of  view 
of  171.99°,  and  azimuthal  angle  of  0°,  and  for  solar  zenith  angles  at 
the  sub-satellite  point  of  0°  and  10°.  It  was  found  that  contributions 
to  the  scattered  intensity  from  collisions  greater  than  20  (that  used 
in  the  Monte  Carlo  calculations)  would  not  contribute  significantly  to 
the  total  scattered  intensity. 

Figure  4 shows  similar  type  data  for  the  case  where  the  cloud 
was  loaded  with  a dust  content  equal  to  1000  times  the  standard  dust 
amount.  The  results  shown  are  for  solar  zenith  angles  0^  - 0°  and  30° 
and  the  receiver  was  looking  in  a direction  defined  by  0 ■ 179.11°  and 
4>  « 0°.  The  scattered  intensity  falls  off  approximately  exponentially 
with  collision  number  and  the  contribution  from  collisions  greater  than 
20  was  estimated  to  be  negligible  for  those  cases. 

The  variation  of  the  reflected  sunlight  intensity  at  the  satellite 
when  the  satellite  zenith  angle  of  view  is  179.11°  is  shown  in  Fig.  5 for 
cases  where  the  azimuthal  direction  of  view,  $,  is  either  0°  or  90°. 

The  data  are  plotted  in  Fig.  5 versus  the  solar  zenith  angle  6q  at  the 
sub-satellite  point.  , It  is  seen  that  there  is  clearly  a difference  in 
the  magnitude  of  the  reflected  radiation  when  the  cloud  is  loaded  with 
100  times  and  1000  times  the  standard  dust.  The  statistics  of  the 
Monte  Carlo  data  for  the  cases  where  the  cloud  contained  no  dust  and 
where  it  contained  the  standard  dust  loading  are  not  good  enough  to  say 
with  certainty  that  a detector  would  see  a difference  in  the  reflected 
radiation  for  those  cases. 
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Fig.  4.  Reflected  Intensity  for  0 - 179.11°  and  <J>  ■ 0 


X STANDARD  DUST 


NO  DUST 


SOLAR  ZENITH  ANGLE  (DEG) 


Fig.  5.  Reflected  Sunlight  for  6 - 179.11 


Figure  6 shows  similar  data  for  reflected  sunlight  when  the 
receiver  zenith  angle  of  view  is  171. 99°.  For  this  receiver  zenith 
angle  of  view,  it  is  seen  that  there  is  a significant  difference  in 
the  magnitude  of  the  backscattered  intensity  for  the  cases  where  the 
cloud  contains  100  and  1000  times  the  standard  dust  loading.  As 
pointed  out  in  the  previous  figure,  it  is  doubtful  that  the  presence 
of  standard  dust  could  be  determined  with  respect  to  the  reflected 
radiation  from  clouds  with  no  dust  or  clouds  loaded  with  standard  dust. 

For  this  satellite  zenith  angle  of  view  and  for  a azimuthal  view  angle 
of  p ■ 0#,  a satellite  would  probably  have  difficulty  determining  if 
the  cloud  contained  as  much  as  100  times  standard  dust  when  the  solar 
zenith  angle  0^  is  greater  than  about  45°. 

A tabulation  of  the  air  emission  signal  at  the  satellite  in 
the  2.653  to  2.760  vim  wavelength  interval  as  a function  of  dust  content 
and  satellite  zenith  angle  of  view  is  shown  in  Table  II.  The  upper 
part  of  Table  II  shows  the  total  intensity  from  all  of  the  atmosphere 
in  the  direction  being  viewed  and  the  lower  part  of  Table  II  shows  the 
intensity  from  air  emission  in  only  the  atmospheric  layer  containing 
the  cloud.  In  general,  the  Infrared  energy  at  the  satellite  in  a given 
viewing  direction  does  not  change  when  either  standard  dust  or  100 
times  standard  dust  is  added  to  the  cloud.  Increasing  the  dust  content 
in  the  cloud  to  1000  times  standard  du9t  does  result  in  a slight  decrease 
in  the  air  emission  reaching  the  detector  in  a given  viewing  direction 
over  that  obtained  for  the  other  dust  contents  considered.  In  general, 
only  about  30Z  of  the  energy  at  the  detector  in  a given  viewing  direc- 
tion comes  from  the  cloud  layer. 

Table  III  shows  that  the  total  intensity  of  air-emitted  radiation 
at  the  detector  in  a given  direction  from  only  the  cloud  layer  is  approxi- 
mately a factor  of  4 to  6.7  larger  than  the  intensity  of  air-emitted 
radiation  which  gets  to  the  receiver  from  the  cloud  without  undergoing 
scattering  in  the  cloud  or  the  upper  atmosphere. 
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Fig.  6.  Reflected  Sunlight  for  0 - 171.99 


TABLE  II.  ENER6Y  AT  SATELLITE  FROM  AIR  EMISSION 
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Standard  Dust  1,941-5  2.043-5  1.913-5  1.844-5  1.827-5 
102  Stand.  Dust  1.893-5  2.073-5  2.083-5  2,007-5  2.011-5 
105  Stand.  Dust  1.250-5  1.512-5  1.552-5  1.699-5  1.699-5 


TABLE  III.  IMPORTANCE  OF  SCATTERING  FOR  AIR  EMISSION 
IN  THE  2,7  un  BAND  FROM  CLOUD  ONLY 

(direct  plus  scattered/direct) 


Dust 

Satellite  Viewing  Angle  (deg) 

Content 

171.99 

173.77 

175.55 

177.33 

179.11 

No  Dust 

6.3 

6.7 

6.5 

6.5 

6.5 

Standard  Dust 

6.1 

6.6 

6.5 

6.5 

6.5 

102  Stand. 

Dust 

5.7 

5.9 

6.1 

5.9 

6.0 

10*  Stand. 

Dust 

3.7 

9.0 

3.9 

4.1 

4.2 
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The  results  obtained  from  the  FLASH  calculations  for  both  the 
reflected  sunlight  and  the  air-emission  radiation  at  the  satellite 
indicate  that  it  is  doubtful  that  one  would  be  able  to  determine  that 
a cloud  has  been  loaded  with  an  amount  equal  to  that  assumed  for 
standard  dust  when  the  detector  is  viewing  infrared  radiation  in  the 
2.653  to  2.760  ym  wavelength  range. 

Since  the  results  shown  in  Figs.  5 and  6 and  Table  II  were 

obtained,  we  have  performed  some  additional  MIE-2  calculations  in  an 

effort  to  determine  what  infrared  radiation  wavelength  interval  offers 

the  best  possibility  of  being  able  to  detect  the  presence  of  dust  in 

high  altitude  clouds.  In  the  preliminary  study  described  above,  the 

cumulus  cloud  particles  were  assumed  to  be  water  droplets.  Since  the 

cloud  was  assumed  to  be  located  between  14  and  15  km  altitude,  where 

the  atmospheric  temperature  is  a minimum,  we  decided  to  assume  that 

the  cloud  was  made  up  of  ice  particles.  The  indices  of  refraction  for 

ice  particles  were  taken  from  Ref.  7.  Figure  7 shows  the  absorption, 

scattering,  and  extinction  cross  sections  computed  with  Mie  theory  for 

the  cumulus  cloud,  assuming  ice  particles.  The  optical  thickness  of 

the  cloud  at  a wavelength  of  0.55  ym  is  16.6  for  a 1.0  km-thick  cloud. 

Figure  8 shows  the  absorption  and  extinction  coefficients  for  the 

standard  amount  of  dust  for  cases  where  the  dust  distribution  is  defined 
-3.5  -4  5 

either  with  a r ’ or  an  r * distribution.  In  both  cases,  the  total 

-11  3 3 

volume  of  dust  was  taken  to  be  5.0X10  cm  /cm  . It  is  noted  in 
Fig.  7 that  a minimum  in  the  extinction  coefficient  for  the  cloud  occurs 
at  10  ym.  It  is  also  to  be  noted  in  Fig.  8 that  a maximum  in  the  absorp- 
tion and  extinction  coefficients  for  Sahara  dust  also  occurs  at  wave- 
lengths near  10  ym.  At  a 10-ym  wavelength  the  absorption  optical  thick- 
ness in  the  cloud  with  standard  dust  will  be  greater  than  it  was  at 
2.7  ym  but  the  total  optical  thickness  will  be  less.  The  emission  from 
the  cloud  with  dust  is  dependent  on  the  magnitude  of  the  absorption 
cross  section  and  the  wavelength.  The  blackbody  emittances  as  given 
by  Planck's  law  for  temperatures  of  207. 5°K  and  276°K  are  shown  in 
Fig.  9 as  a function  of  wavelength.  The  207. 5°K  temperature  is  the 
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average  temperature  of  the  cloud  and  276°K  Is  the  maximum  temperature 
which  occurs  at  50  km  altitude.  For  the  2.7  |im  wavelength,  Planck's 
law  shows  that  there  is  a factor  580  more  energy  emitted  at  a tempera- 
ture of  276°K  than  there  is  at  a temperature  of  207. 5°K.  Increasing 
the  wavelength  of  emission  to  10.0  Urn  results  in  a much  larger  increase 
in  the  magnitude  of  the  blackbody  emittance  as  given  by  Planck's  law. 

In  addition,  the  emittance  at  a wavelength  of  10.0  pm  for  a temperature 
of  276°K  is  only  a factor  of  5.88  greater  than  the  emission  at  a wave- 
length of  10.0  pm  when  the  atmosphere  is  at  a temperature  of  207. 5°K. 

In  the  next  few  months,  RRA  plans  to  run  FLASH  calculations  at 
other  wavelengths  for  both  reflected  solar  radiation  and  air  emission 
to  see  if  there  are  wavelengths  for  which  it  is  possible  to  differentiate 
between  clouds  that  contain  either  standard  dust  or  no  dust. 


21 


REFERENCES 


Wolfram  G.  Blattner,  Utilization  Instructions  for  Operation  of 
the  MIE-2  Program  on  the  Sandla  CDC-6600  Computer,  Radiation 
Research  Associates,  Inc.  Research  Note  RRA-N7513  (September 
1975) 

Wolfram  G.  Blattner,  Utilization  Instructions  for  the  FLASH 
Monte  Carlo  Frocedure,  Radiation  Research  Associates,  Inc. 
Research  Note  RRA-N7512  (August  1975) 

R.  A.  McClatchey,  private  communication,  1974 

F.  Volz,  "IR  Optical  Constants  of  Ammonium  Sulphate,  Sahara 
Dust,  Volcanic  Pumice  and  Flyash",  Appl.  Opt.,  Vol.  12,  No.  3, 
p.  564  (1973) 

D.  Deirmendj ian,  Electromagnetic  Scattering  on  Spherical  Poly- 
dispersions , Elsevier,  New  York  (1969) 

G.  M.  Hale  and  M.  R.  Querry,  "Optical  Constants  of  Water  in 
the  200  nm  to  200  pm  Wavelength  Region",  Appl.  Opt.,  Vol.  12, 
No.  3,  p.  555  (1973) 

W,  Irvine  and  J.  Pollack,  "Infrared  Optical  Properties  of  Water 
and  Ice  Spheres",  Icarus  (1967) 


No  of 
Copies 

2 

2 

12 

2 

5 

1 

3 

1 


DISTRIBUTION 


No  of 
Cop  ios 


DOD  AGENCIES  1 

Di rector 

Defense  Nuclear  Agency 
ATTN:  RAAL  1 

Technical  Library 

Di rector 

Defense  Intelligence  Agency 
ATTN:  Nuclear  Energy  Division 
Technical  Library 

1 

Defense  Documentation  Center 
ATTN:  Document  Control 

2 

Di rector 

Defense  Advanced  REsearch 
Projects  Agency 

ATTN:  NMRO  1 

TIO 


ENERGY  RESEARCH  AND 

DEVELOPMENT  ADMINISTRATION  2 

Los  Alamos  Scientific  Laboratory 

ATTN:  J-10  (Dr  H.  Horak) 

J-10  (Dr  Guy  Barasch) 

T-12  (Dr  D.C.  Cartwright) 

TD-3  (Dr  Bob  Henson)  3 

Technical  Library 

ERDA/DMA 

ATTN:  Mr  D.  Gale 

Sandia  Laboratories  1 

ATTN:  Dean  Thornbrough 
Bob  Bradley 
Reports  Collection 

Lawrence  Livermore  Laboratory  2 

ATTN:  Technical  Library 

Oak  Ridge  National  Laboratory 

ATTN:  Technical  Library  23 

QOt*iGD08s0(?3@(s] 


EG&G,  Inc. 

Bedford  MA 

ATTN:  Technical  Library 

EG&G , I nc . 

Los  Alamos  NM 
ATTN:  Dr  C.  Mitchell 


DEPARTMENT  OF  T I IE  ARMY 

Army  Nuclear  Effects  Lab 
ATTN:  Technical  Library 

US  Army  Electronics  Command 
ATTN:  Technical  Library 
Dr.  R.  Buser 

US  Army  Combat  Developments 
Command-Insti tute  of  Nuclear 
Studies 

ATTN:  Technical  Library 

US  Army  Waterways  Experiment 
Station-Mobility  and 
Environmental  System  Lab 
ATTN:  Dr.  L.E.  Link 

Dr.  Warren  Grabau 

US  Army  Electronics  Command 
Atmospheric  Sciences  Laboratory 
ATTN:  Dr  M.L . Vatsia 
Dr  Jerry  Lentz 
Dr  Richard  Gomez 

US  Army  Missile  and 
Development  Command-Physical 
Sciences  Directorate 
ATTN:  Dr  Julius  Lilly 

Ballistic  Research  Laboratory 
ATTN:  Mr.  Norman  Banks 
Mr.  Ennis  Quigley 


1 


No  of 
Copies 

1 

1 

2 

1 

4 

1 

4 

1 

2 

1 


®fflGDas§8(?9©d] 


No  of 
Copies 


DEPARTMENT  OF  THE  NAVY 

Naval  Research  Laboratory 
ATTN:  Technical  Library 

Naval  Postgraduate  School 
ATTN:  Technical  Reports 
Librarian 

Naval  Surface  Weapons  Center 
ATTN:  Mr  Williard  Derksen 
Mr  Barry  Katz 

Naval  Ordnance  Test  Station 
ATTN:  Dr.  Richard  Bird 


DEPARTMENT  OF  THE  AIR  FORCE 

AF  Geophysics  Laboratories 

ATTN:  OPR  (Mr  H.  Gardiner) 
OPA  (Dr  R.  Fenn) 

OPR  (Dr  A.  T.  Stair) 
Research  Library 

AF  Institute  of  Technology 

ATTN:  Library 

AF  Weapons  Laboratory 

ATTN:  Dr  Joe  Janni 
Mr  A1  Sharp 
Dr  C.  Needham 
Technical  Library 

Space  and  Missile  Systems 
Organization 

ATTN:  SZS  (Maj  H.  Hayden) 

AF  Technical  Applications 
Center 

ATTN:  TFR  (Capt  J.  Lange) 


AFETAC 

ATTN:  DNE  (Capt  L.  Mendenhall) 


POD  CONTRACTORS 

3 Aerospace  Corporation 

ATTN:  Dr  R.D.  Rawcliffe 
Dr  Kelly  Spearman 
Library  Services 

1 Analytical  Systems  Corporation 

ATTN:  Dr  I.  Kohl  berg 

1 Battel le  Memorial  Institute 

ATTN:  Radiation  Effects 

Information  Center 

1 AVCO  - Everett  Research 

Laboratory 

ATTN:  Technical  Library 

1 General  Electric  Company-TEMPO 

ATTN:  DNA  Information  and 
Analysis  Center 

1 Kaman  Sciences  Corporation 

Kaman  Nuclear  Division 
ATTN:  Technical  Library 

1 Lockheed  Missiles  and  Space 

Company 

ATTN:  Technical  Library 

1 The  MITRE  Corporation 

ATTN:  Technical  Library 

1 The  RAND  Corporation 

ATTN:  Technical  Library 

1 R&D  Associates 

ATTN:  Technical  Library 

1 Stanford  Research  Institute 

ATTN:  Technical  Library 


24 

QDod<3D<2iss8(?0©(£I 


f 


(DffleDassfitfted] 


No  of 
Copies 


1 

1 

1 

1 

1 

1 

20 

1 

1 


TRW  Systems  Group 

ATTN:  Technical  Information  Center 

Visidyne,  Inc. 

ATTN:  Dr.  J.  Carpenter 

Mission  Research  Corporation 
ATTN:  Technical  Library 

Teledyne  Brown  Engineering 

Science  Applications,  Inc. 

ATTN:  Dr  R.  Hillendahl 

HSS,  Inc. 

ATTN:  Dr  H.  Stewart 

Radiation  Research  Associates 
ATTN:  M.B.  Wells 

Bendix  Engineering  Corporation 
Grand  Junction  Operations 
ATTN:  Dr  Peter  Alexander 

Institute  for  Defense  Analyses 
ATTN:  Technical  Information  Office 


25 


